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Solitons confined in channels are studied in the two-dimensional nonlinear Schrodinger equation. 
We study the dynamics of two channel-guided solitons near the junction where two channels are 
merged. The two solitons merge into one soliton, when there is no phase shift. If a phase difference is 
given to the two solitons, the Josephson oscillation is induced. The Josephson oscillation is amplified 
near the junction. The two solitons are reflected when the initial velocity is below a critical value. 
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' The one-dimensional nonlinear Schrodinger equation has been intensively studied as a typical soliton equation for 
,. optical solitons in optical fibers [1] and matter- wave solitons in the Bose-Einstein condensates (BECs) [2], [3]. However, 
solitons in two- or three-dimensional nonlinear Schrodinger equations were not intensively studied. Two-dimensional 
solitons exist in guided channels in the two-dimensional nonlinear Schrodinger equation. [4] The guiding channel can 
^rfH ' be constructed by modifying the profile of the refraction index in an optical planar waveguide. Cigar-shaped traps are 
^ ] used as a guiding channel to confine matter-wave solitons in BECs. Solitons can propagate along the guiding channel 
■ with an arbitrary velocity if the guiding channel is uniform and the norm of the sohtons is below a critical value for 
' collapse. In our previous study [5], we investigated the reflection of a channel-guided soliton in a tapered channel and 
the splitting at a branching point where one channel branches into two channels [6]. In this study, we investigate the 
p ^ motion of solitons in curved channels and near the junction where two channels merge into a channel. 
The model equation is written as 
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i^ = -\v^cj,-\<j,\^<j, + U{x,y)(t>, (1) 



where U{x, y) denotes the potential for confinement, which is used to make a guiding channel. For a straight channel, 
^ , U{x,y) — —Uq for —xq < x < xq and U{x,y) = for other regions. The width of the channel is denoted as 2xq and 
■ the depth of the potential is denoted as Uq. There are stationary soliton solutions 4)o{x,y) in the straight channel. 
CO I The stationary solution 4>o(x, y) can be numerically obtained by the imaginary-time evolution of eq. (1) as was carried 
\l I out in our previous study 0. The two-dimensional soliton can move with any velocity ky if the initial condition is 
set to be 4>Q{x,y - yo)exp{iky{y - ya)). 
, In the previous study, we investigated the splitting and reflection of a two-dimensional soliton in a branching 
' channel and designed a branching channel where a soliton can split smoothly into two solitons without changing its 
, velocity. If the branching channels merge into a straight channel at a junction point again, a waveguide system, such 
as the Mach-Zehnder interferometer, can be constructed, as shown in Fig. 1(a). The branching point \s y = —40 
and the junction point is y = 40 in the channel system in Fig. 1(a). A phase shift A(/5 is given to the two split 
solitons when the center of the solitons goes through the line y = as 4'{x,y) — ^ (p{x,y) exp{iALp) for a; < and 
<j){x,y) — > ^(x, y) exp(— iA(/3) for a; > 0. This procedure is an artificial one, which is simple and easy for numerical 
. simulation. In an actual experiment, the phase shift can be effectively given to the two solitons by controlling the 
profiles of the refraction index and the forms of the two channels similar to the typical Mach-Zehnder interferometer. 
When Aiy9 ~ 0, the two solitons merge again into one soliton at the junction located at y = 40. We have performed 
numerical simulation with the split-step Fourier methods with 128 x 1024 modes for this rectangular system. For 
nonzero Atp, a zigzag oscillation of the center of gravity is observed, as shown in Fig. 1(b). Figure 1(b) shows the 
trajectories of the center of gravity of the solitons for A(^ = 0, 7r/64, 7r/32 and 7r/16 at ky = 1. The norm N of the 
soliton is 5. A zigzag oscillation appears in a region of y > 20. The amplitude of the zigzag oscillation increases with 
Aiy9, but the period of the oscillation is almost the same. This zigzag oscillation is interpreted as a remnant of the 
Josephson oscillation that appears before the junction point at y = 40. The Josephson oscillation occurs owing to the 
tunnel effect through the potential wall of height Uq between the two channels. The Josephson oscillation is explained 
in detail in the following. 

When the initial velocity ky is decreased for a fixed phase shift value, the solitons are reflected before the junction 
point at y — 40. Figure 1(c) shows the time evolution of |(/)| along the center lines of the left and right channels at 
ky = 0.37 for A(p = 7r/256. Before the reflection, the amplitudes of the two solitons in the left and right channels 
are almost the same, but the difference between the amplitudes of the two solitons increases after the reflection, and 
the soliton becomes localized in one of the two channels. The critical velocity kyc is 0.38 for Aip = 7r/256. The 
relationship between kyc and the phase shift A(p is shown in Fig. 1(d). Note that the horizontal axis is plotted with 



2 



(a) (b) 




FIG. 1: (a) Channel with branching and junction points, (b) Trajectories of the center of gravity of sohtons for ky = 1. The 
phase shifts of Aip = 0, 7r/64, 7r/32 and n/16 are given at y = 0. The straight hne is the trajectory for Aip = 0. The zigzag 
oscillation grows with Aip. The trajectory that exhibits the largest zigzag oscillation corresponds to Aip = 7r/16. (c) Time 
evolution of 101 along the center lines of left and right channels at ky — 0.37 for Aip = n/256. (d) Critical values kyc of the 
reflection as a function of A^p. 
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FIG. 2: (a) Time evolution of \u\ (solid curve) and \v\ (dashed curve) for ky = 0.29 and A^p = 7r/256. (b) Time evolution of 
R = {Nu - Ny)/{Nu + iV„). (c) Time evolution of i{t) determined using eq. (7) for ky = 0.29 and AO = 7r/128. (d) Time 
evolution of R determined using eq. (7). 



a logarithmic scale. Figure 1(d) implies that the critical velocity kyc increases rapidly near Atp — 0. 

If there is a phase difference between two sohtons located on different channels, the Josephson oscillation oc- 
curs [7], [8], because the two solitons interact with each other by the tunnel effect through the potential barrier of height, 
Uo- We can take an ansatz for the form of as (f>{x, y) = u{y, t) exp[— {a; — ?7(y)}^/(2&^)] exp{—ifit) + v{y, t) exp[— 
77(?/)}^/(26^)] exp(— i/xt), where i'ri{y) is the x-coordinate of the central point in the right and left channels at 
y. The substitution of the form into the Lagrangian (2) and the variational principle: d/dt{5L/Sut) = 5L/5u, 
d/dt{dL/Svt) — dL/6v, yield approximately coupled equations for u and v: 

2dy^^ c(|u|2 + g\v\^)u - d{v - u), 

l^-ci\v\^+g\u\^)u-d{u-v), (2) 

where g = 2exp{—2'ri{y)^ /b^) and d — 2Uo exp{—2r]{y)^ /b'^). Here, we have neglected some complex terms, such 
as u*w^ and The parameter c is assumed to be c = l/{\/2(l + g)}, because the symmetric solution u — v 

propagates without changing the profile of the solitons. When rj is large, mutual interaction is weak and the two 
solitons propagate independently. However, as y is close to the junction point y = 40, the Josephson effect becomes 
strong because of a small 77 value. 

Figure 2(a) shows the time evolution of |m| and \v\ determined using the coupled equations (2) at ky — 0.2. 
The initial conditions of u and v are u = ^/cosh(?;/W) exp(zA(^) and v = A/cosh{y/W)exp{—iAip) where b = 
0.9, A^o = 5/{^/nb),W = AV2/N0, A = ^/No/{4W), and Aip = 7r/256. The two solitons are reflected, and one of the 
solitons dominates after the reflection. The critical velocity ky = 0.295 is slightly different from that in the direct 
numerical simulation, but the qualitative behaviors are similar. The asymmetry of the two solitons is expressed as 
R = {N^ — Ny)/{Nu + Ny) where — J \ u\'^dy and Ny = J \v\'^dy are the norms of u and v, respectively. The time 
evolution of the ratio R is shown in Fig. 2(b). The ratio R increases from and reaches a positive constant value 
after the reflection. This implies that one soliton dominates in the u-channel after the reflection. 
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liu{y,t) = Asech{{y — £_) /W} exp{ip{y — — i9i} a.ndv{y,t) = Bsech{{y — ^) /W} exp{ip{y — ^) — i62} are further 
assumed, the effective Lagrangian is evaluated as 

^eff ~ \ J {^(^t^* ~ ""t + ""tv* - vVf) - \uy\'^ - \vyf + c(|u|^ + + 2fi(|M|^|t;p) 

+2d{uv* + vu*) - 2d(|up + \vf)}dy 
= NoHt - pV2 - l/i&W^) + {Oit + ^2t)/2 + RiOu - ^2t)/2 

+R'^Noe/{16W'^) + rf{ Vl - cos(6i2 - Oi) - 1}] (3) 

where No = J{\u\^ + \v\'}dy = 2{A^ + B^)W , R={A^- B^)I{A^ + B% e(0 = / c{y){l - g{y)}sec\i\{y - ^)/W}dy, 
and d{£) = J d{y)sech^{{y — ^)/W}dy/ J sech^{{y — ^)/W}dy. The variational principle yields 
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w = 



{1 + R^ + {1- R^)g}Noc' 
P, 



dt 

dp R'^Node , r ,^ ^ . ^-.dd 

di = i6H^a! + {^^^°^('^-'^)-^>af' 

dR 
dt 

dAO RNoe 2dR 



2JVl-i?2sinA6l, 

cosA6i, (4) 



dt VT^^ 

where A6 — 02 — 01, and ^ is the y-coordinate of the center of gravity of the solitons. If i? = and di = O2, 
dp/dt — and the two solitons propagate with a constant velocity. If R << 1 and AO << 1, d^A9/dt^ — —2d{2d — 
Noe/{4W'^)}A9 is obtained using the last two equations in cq. (4), which describes the Josephson oscillation. If 
Nqc / {AW'^) — 2(i > 0, the Josephson oscillation is amplified and the symmetric state u = v becomes imstablc. Because 
e(^) and (i(^) arc not uniform, dp/dt can become negative owing to the amplified Josephson oscillation. Figures 2(c) 
and 2(d) show the time evolution of ^ and R for A^o = 5/{^/TTb) with b = 0.9 and ^(0) = 0,p{0) = ky = 0.2, and the 
initial phase difference A9{0) = 2Aip = 7r/128. The reflection of the trajectory and the amplification of the Josephson 
oscillation are observed. The ratio R takes a constant value for large t values, although the sign of the constant 
value depends strongly on the initial velocity p(0) because of the many oscillations of R near the reflection point. 
The critical value kyc for A^(0) = 7r/128 determined using eq. (4) is 0.305, which is close to that by eq. (1). The 
variational approximation is fairly good. 

To summarize, we have performed some numerical simulations of channel-guided solitons in the two-dimensional 
nonlinear Schrodinger equation. We have found that reflection occurs near the junction induced by the amplification 
of the Josephson oscillation. The complex dynamical behaviors can be analyzed approximately on the basis of the 
variational principle using effective Lagrangians. Channel-guided solitons exhibit various complex dynamics, and we 
would like to investigate this further in the future. 
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